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A TECHNIQUE FOR PEAK POWER TO AVERAGE POWER RATIO 
REDUCTION IN MULTICARRIER COMMUNICATION SYSTEMS 



Inventors 

Amit Verma and M T Arvind 

FIELD OF INVENTION 

This invention pertains to communication systems that employ a multitude of carriers for 
conveying information through a transmission medium. In particular, this invention addresses the 
problem of peak to average to power ratio reduction in such systems. 

BACKGROUND OF THE INVENTION 

A multicarrier transmission system is one that employs Frequency Division Multiplexed 
(FDM) subcarriers for transmission of data across a communication channel. 

A typical multicarrier transmission system is depicted in Figure 1. It shows serial data in 
the form of bits input to a serial to parallel converter 100, which also frames the bits. (A frame is 
an ordered sequence of bits of a given size). The output bit frame is mapped to symbols in the 
mapper 300. The frame of symbols then passes through the modulator 400. The modulation 
could be carried out, for instance, using Inverse Fast Fourier Transform (IFFT) and the 
corresponding demodulation carried out using Fast Fourier Transform (FFT). The modulated 
signal is converted to serial data (block 700) and transmitted on to the channel 900, where it 
suffers several impairments, in addition to being corrupted by noise. At the receiver, the received 
data is grouped into frames by a serial to parallel converter 1000 and subsequently passed 
through the demodulator 1100. The demodulated symbols are converted to bits through an 
inverse symbol mapper 1300, which also performs parallel to serial conversion. 
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The data bits input to the transmitter module are independent and the symbols obtained 
by the bits-to-symbol mapper are also independent. Whenever the number of subchannels is 
large, the modulator output magnitude tends to have a Gaussian distribution, which has an 
infinite peak to average power ratio. A clipper 800, which has a cutoff based on a predesignated 
5 clipping probability, is used to limit the peak signal power that is transmitted on to the channel. 
The clipping probability is kept low so that the error introduced due to clipping is small. But if 
the predesignated clipping probability is reduced, the increase in peak power may cause 
saturation in the subsequent amplification stages and intermodulation distortion due to 
nonlinearities in the transmission medium. Thus it becomes necessary to limit the peak power 

V30 transmitted on to the channel by reducing the clipping level while still maintaining a low 

^ clipping probability. 

sj Figure 1 shows a peak detector 500 at the output of the modulator, which activates a 

symbol modifier 600 whenever the magnitude of any of the modulator output samples exceeds a 
predesignated threshold or set of thresholds. Although Figure 1 shows the symbol modifier 

;g5 correcting the symbols that are being modulated, it is possible to achieve the same effect by 

3 

correcting the modulated samples. The threshold could be set based on the level of peak power 
that can be transmitted on to the channel, while the task of the symbol modifier is to modify the 
symbols being modulated in such a manner that the peak is reduced. The symbol modifier needs 
to take care of a number of objectives: 
20 1. The clipping probability should remain unaltered when the clipping threshold is 

lowered. 

2. The loss in data rate caused due to symbol modification should be minimal. 
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3. A preferable requirement is that the symbol modifier should have a low complexity. 

4. The PAR reduction scheme should be transparent to the receiver. 

5. The modifications made to the data symbols should not cause an increase in the 
average subchannel powers. 

6. The modifications introduced in the symbols should not cause significant increase in 
the error rate at the receiver end. 

Most of the above requirements are conflicting and it may not be possible to meet all of 
them at the same time. The amount of importance given to each of these requirements by a 
particular method decides the amount of Peak to Average power Ratio (PAR) reduction 
achievable by that method. 

Most of the prior art addresses several of the requirements listed above, but not all of 
them together. Some of the conventional methods use Peak Reduction Tones (PRT) (a 
predesignated set of subchannels) to reduce the PAR. One of the PRT methods is described by J. 
Tellado and J. Cioffi in "PAR Reduction in Multicarrier Transmission Systems (97-367)", 
T1E1. 4/97-367 December 8, 1997. It uses an iterative procedure which aims at reducing the 
largest peak every iteration. Another method is described by the same authors in "Revisiting 
DMT's PAR (98-083)", T1E1. 4/98-083, March 3, 1998. It uses an iterative procedure that tries 
to reduce all the peaks present during each iteration. Both the methods employ a precomputed 
peak reduction kernel. The reduction in PAR is obtained at the loss of several data carrying 
subchannels. PRT methods require that the receiver be aware of the predesignated subchannels to 
be used for peak reduction. 

The modulo-D scheme described in "A new approach to PAR control in DMT Systems", 
NF-83, ITU-T, Study Group 15, 11-14 May 1998 uses an expanded constellation in each 



subchannel for PAR reduction. The method uses symbols drawn from larger constellations 
whenever the original modulation produces large peaks. The use of larger constellations results 
in increased average power and higher intermodulation distortion. The receiver will have to 
perform a reverse operation to the transmitter action to recover the original data symbols. 
5 A sign inversion method has been specified in "An Efficient Implementation of PAR 

reduction method based on subset inversion", AB-061r2, ITU-T, Study Group 15, 3-7August, 
1998. It consists of dividing the set of subchannels into several subsets and achieving a PAR 
reduction by inverting the sign of the symbols in one or more of the subsets. The transmitter uses 
part of the data carrying bits to indicate to the receiver whether a sign inversion has been applied, 
fl 0 This method requires receiver cooperation and results in data rate loss. Additionally, any error in 
; ! *f a subset inversion indicator bit will result in incorrect decoding of data bits associated with that 
Q subset. 

SUMMARY OF THE INVENTION 

This invention proposes a peak-sensitivity based scheme for introducing modification(s) 
to the frame of data symbols, such that the PAR is reduced. The method employs a penalty 

. 

function that assumes a non-zero value only when there is at least one sample whose magnitude 
exceeds a predetermined peak threshold. Whenever the penalty function becomes non-zero, a 
sensitivity-based scheme is proposed for computing the corrections to the transmit symbols. The 
amount of correction provided by the method is flexible. It is based on the tradeoff between the 
20 PAR reduction required and the error margin available for symbol decoding in individual 
subchannels at the receiver. The advantage of the proposed scheme is that the receiver need not 
have prior knowledge of the transmitter employing a PAR reduction mechanism. The receiver 
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need not perform any extra action that reverses the effect of the transmitter action for achieving 
PAR reduction. No channels need to be specifically allocated for PAR reduction signals. 

The advantage of the scheme proposed is that it satisfies requirements 1, 2, 3, 4 and 5 
completely and allows for tradeoff in requirement 6 (in terms of compromise between PAR 
reduction and increase in symbol errors at the receiver). Several modifications to the proposed 
scheme are also suggested in order to reduce the computation overheads. The schemes involve 
usage of a precomputed peak kernel, which is scaled and rotated suitably such that the peaks are 
reduced. The main difference between the peak kernel proposed in this invention and the one 
used in prior art is that the proposed peak kernel could have non-zero values in all the 
subchannels. But the individual values in each subchannel are chosen to be small, so that 
negligible deterioration in SNR takes place at the receiver and there is little or no increase in 
error rate whenever peak correction is applied. 

The salient features of the PAR reduction scheme proposed in this invention are: 

1 . The scheme modifies the transmit data frame to obtain PAR reduction, while having a 
small effect on the probability of error. The modifications could be performed on 
either the unmodulated or the modulated samples. 

2. The scheme has low complexity and can be implemented efficiently. 

3. The amount of PAR reduction can be traded off against increase in error rate. 

4. The transmitter actions with respect to PAR reduction are totally transparent to the 
receiver, since the receiver actions for recovering the transmitted data remain 
unaltered by this scheme. 

5. Modifications to the proposed scheme are also described, which reduce the 
computation while maintaining a good PAR reduction. 



BRIEF DESCRIPTION OF THE DRAWINGS 



Figure 1 is a block diagram of a typical multicarrier system that employs a forward error 
correction code. 

Figure 2 is a block diagram of the proposed scheme. 

Figure 3 is a flow diagram of the gradient sensitivity based PAR reduction method. 
Figure 4 depicts an iteration of Procedure 1. 
Figure 5 describes Procedure 2. 

DETAILED DESCRIPTION OF THE INVENTION 

The proposed scheme is applicable to multicarrier communication systems where the 
input data bits are converted to a symbol modulated by each subcarrier, modulated and 
transmitted on the channel. More generally, it is applicable to any system where the modulation 
and demodulation are performed at the frame (symbol) leveL However, the ensuing description 
considers a multicarrier system employing IFFT and FFT for modulation and demodulation, 
respectively. 

A block schematic of the proposed method is shown in Figure 2. The data frame of M 
bytes is mapped to symbols in block 2200 and modulated using the IFFT block 2300. It is 
assumed that the multicarrier system has N/2 channels. The output rate is assumed to be N 
samples per frame (i.e., the number of IFFT samples per frame is N). 

The IFFT samples are passed through a peak detector 2400. The peak detection is carried 
out on the modulated samples using a predetermined threshold (or set of thresholds). The output 
of the peak detector activates the symbol modifier block 2500 whose output is added to the 



modulator output. This is finally passed through the Digital to Analog converter, filters and line 
drivers (represented in block 2700) and transmitted onto the channel. 

The functions of peak detector and symbol modifier blocks are explained below. 
The peak detector 

The inputs to the peak detector are N samples of the IFFT. The peak detector uses a 
predesignated threshold T. Each IFFT sample is compared with the threshold and the indices of 
samples whose magnitude exceeds the threshold are stored along with the amount by which they 
violate the threshold. This forms the input to the symbol modifier block. 
The symbol modifier 

This block determines the modification(s) to be made to the input symbols (or samples) 
such that the output of the symbol modifier reduces the PAR of the encoded frame. Below 
several schemes are disclosed of varying degrees of complexity. The clipping level is fixed based 
on the permissible error rate introduced by the transmitter module. In general, the clipping 
probability for a given threshold is small. In such a case, the probability of more than one sample 
being clipped in a frame is small compared to the probability of a single sample (peak) getting 
clipped. 

The gradient sensitivity based approach for computing the symbol modifier is now 
described (and depicted in Figure 3). In this method, a differentiable penalty function is used to 
compute a measure of peak threshold violations (block 3000). Partial derivatives of the function 
with respect to subchannel symbol magnitudes and phase angles is computed to determine the 
sensitivity of the peaks with respect to the data symbols (block 3100). The gradient vector is used 
to compute the direction of variation of the symbol magnitudes and phases (block 3200). The 
optimal amount of variation in each direction can be computed using traditional optimization 



techniques, but the computation overhead for this operation can be very high. Further, the 
modifications suggested for the best possible PAR reduction can cause a large increase in the 
error rate, offsetting the gains obtained by performing PAR reduction. A compromise solution 
results by considering only the direction given by the gradient vector and making small 
modifications to a few symbols, such that the PAR is reduced (block 3300). 

Denote the frame of data symbols by X = (r 0 , r t exp(J9i) ? r 2 exp(j0 2 ), r W2 _ x exp(jO N/2 _ 
i)? r N/2)> where r { and Q { denote the magnitude and phase of symbol in channel i. Note that the dc 
value and Nyquist channels have real values. The IFFT operation yields time domain samples 
x[k],k-0,l, ...,N-1; where 



N/2-1 



x[k] = r 0 + 2 ^ r- cos 



2nki 

If 



+ r N/2 cos(ftk) 



The outputs of the peak detector are the locations where the IFFT sample magnitudes 
violate the threshold T and the extent of the violation. Thus, for a sample x[k] > T, the violation 
is indicated as (x[k]-7), while it is (x[k]+7) if x[k]+r<0. The following penalty function is used 
for each point x[k] that violates the threshold T. 



h(x[k]) = 



(x[k]-T) 2m if x[k]>T 
0 if \x[k]\ < T 

(x[k] + T) 2m if x[k]<-T 



where m is a positive integer that decides the severity of penalty. In the preferred embodiment, 
we have taken m to be unity in order to minimize computation. Denoting the vector (x[0], x[l], 
. . . , x[N-l]) by x, the net penalty function f(x) is computed as the sum of the individual penalty 



N-l 



functions, i.e., f(x) = ^ h(x[k]) 



k=0 
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Although the summation is indicated over all possible indices, the sum will comprise 

only of the non-zero values corresponding to each peak violation. The proposed sensitivity based 
method involves computation of partial derivatives of f(x) with respect to each magnitude and 
phase angle. 



Note that when m=l, the derivative of h(.) corresponds to the amount by which a peak 
threshold is violated. The objective of locally minimizing the multi-modal penalty function f(.) is 
computationally expensive. Global minimization of f is not feasible since the resulting symbol 
modification may correspond to excessive amount of noise being added to the symbol In order 
to obtain a good PAR reduction with a small amount of computation, the following strategy is 
employed. The negative direction of the gradient corresponds to reduction of PAR. Also note that 
even a small variation in a parameter for which the corresponding term in the gradient vector has 
large magnitude, causes a larger reduction in the peak value. Further, if the requirement of not 
altering the frame power has to be satisfied, only the phase angles can be changed. The amount 
of variation in symbol phases and/or magnitudes should be such that it causes minimal 
deterioration in the SNR of the received symbols. Keeping this in mind, a small fraction (5) of 
each variable could be kept as the limit of change. The first procedure involves changing the 
symbol values by small values to effect a peak reduction. 
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df 




+ {l,...,AT/2-l} 



cos(Vzfc) 



The pseudo code for the procedure is given below: 
Procedure 1 

1. Determine the values 8 0 , S 1? .., 8 N/2 by which the symbols in each subchannel can be 
changed. (The change may have amplitude and phase components.) 

2. Compute the gradient vector as mentioned above. 

3. Apply a correction to the input vector in the opposite direction of the gradient. The 
amount of correction is limited to the value of 5^ for each subchannel. 

4. Repeat steps 2-3 if required. 

5. Remodulate the changed symbols and transmit the modified samples. 

Procedure 1 is depicted in Figure 4. The values of perturbations for each symbol are 
determined (block 4000). The correction to the symbols is computed by using the direction of the 
gradient vector (block 4100). The symbols are modified by this correction, modulated and 
transmitted on to the channel (block 4200). 

The values for 5 Z depend upon the bit to symbol mapping scheme. For example, consider 
a 4-point constellation with points (1,1), (1,-1), (-1,-1) and (-1,1). The inter-point distance is 2 on 
each axis. Hence the maximum value of h { is 1 on each axis so that the modified symbol is not 
mistaken for some other point. In practice, b { is chosen as a fraction of the interpoint distance so 
that the SNR in the particular subchannel is not degraded too much. For the value of 8 Z chosen, 
the decrease in the SNR in the subchannel i at the receiver side can be computed based on the 
clipping probability and the average constellation power. It is clear that the deterioration in the 
SNR in the ith subchannel at the receiver is more when corrections are applied more frequently 
(i.e., T is small and the clipping probability is large) or if |8J is large. Thus the values of |8J and T 
can be chosen based on the tradeoffs between PAR reduction and SNR deterioration in each 
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subchannel on the receiver side. If the gradient with respect to only the phase angles is computed 
and used in procedure 1, the power of the transmit symbol frame remains unaltered. 

If the correction magnitudes are chosen to be large, it is possible that new peaks might 
result while correcting for the existing peaks. In such cases, the corrections can be repeatedly 
applied as indicated in Procedure 1. However, iterative correction needs repeated computation of 
the gradient vector and additional care to avoid repeated and excessive correction in the same 
subchannel. The gradient computation is also computationally expensive. An alternative 
procedure, which does not need gradient computation and uses a pre-computed peak reduction 
kernel to compensate for a single peak, is now described. 
Procedure 2 

1 . Determine the values b { as in procedure 1 . 

2. Compute the peak reduction kernel as the IFFT of vector of b { values in each 
subchannel. 

3. Rotate the peak reduction kernel by an amount such that its peak coincides with the 
peak of the output frame. Multiply the peak reduction kernel by -1 if its peak has the 
same sign as the frame peak. 

4. Add the rotated, possibly sign-inverted peak reduction kernel to the output frame 
samples. 

5. Transmit the modified samples. 

Procedure 2 is depicted in Figure 5. The perturbation values are fixed and the peak 
reduction kernel computed (block 5000). Whenever a peak is detected, the kernel is rotated to the 
location of the peak and added with proper sign (block 5100). The modified data is transmitted 
on to the channel (block 5200). 
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The peak reduction kernel is computed as follows. 

The 5 { values consist of an amplitude and a phase component. The phase component of 
each Sj may be set to zero. This ensures that the peak reduction kernel has its peak value at the 
first sample and that the magnitude of the peak is the maximum attainable using the chosen set of 
5 { values. The rotation of the peak reduction kernel to coincide with a chosen peak of the data 
frame corresponds to rotation of the input 8, sample by an angle proportional to i. 

Multiple peaks are handled in the following manner. For each IFFT sample that violates 
the threshold T ? the peak reduction kernel is rotated as described in procedure 2 to coincide with 
that peak location. Each rotated kernel is scaled by a factor which is a function of the extent of 
peak violation. The scaling factors are chosen such that the sum of their magnitudes is 1. This 
ensures that the noise power thus added is limited. 

The preferred embodiment involves use of an IFFT for modulation. Any alternate means 
of frame (symbol) based QAM modulation may also be used. This could also include PAM and 
DPSK modulation. 

The preferred embodiment uses symmetric, 2-dimensional constellations for mapping bits 
to symbols. The procedure also applies to any subset of the 2-dimensional constellation or multi- 
dimensional constellations. 

The symbol modification can be applied to its phase or amplitude or both. 

Several of the modules can be implemented in hardware or software. 

The preferred embodiment used a specific penalty function that is the summation of even 
powers of each peak violation. Instead of simple summation, suitably weighted summation can 
be used. More generally, the penalty function itself may be different. 
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The gradient directions for single peak occurrences can be pre-computed and stored. The 
computation requirement in this case reduces to one of finding a weighted sum of the directions 
for each peak. 

Several other heuristics known to those of ordinary skill in the art can be used to compute 
the peak reduction kernel of procedure 2 and are considered to be within the scope of the present 
invention. 
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Claims 

1 

2 1 . In a multicarrier communication system in which a signal to be transmitted comprises 

3 data bits to be converted into a symbol modulated by each subcarrier of the signal prior to 

4 transmission on a channel, a method for minimizing a peak to average power ratio while 

5 minimizing introduction of errors into the signal to be transmitted: 

6 sampling the symbols to be transmitted of a frame; 

7 compare magnitudes of the samples of the frame to a predetermined threshold to 

8 determine whether sample magnitudes in the frame violate the 

9 predetermined threshold, the predetermined threshold being selectable to 
flO control the number of samples violating the threshold; 

ps 

H 1 responsive to determining a sample magnitude does violate the predetermined 

*32 threshold, applying a differentiable penalty function to the samples having 

13 magnitudes exceeding the predetermined threshold; 

Iff 4 computing a net penalty function value, the net penalty function value responsive 

; 3 "5 to the individual penalty function values computed for the samples having 

1 6 magnitudes exceeding the predetermined threshold; 

1 7 computing a gradient vector responsive to the net penalty function value; 

1 8 determining a direction of the gradient vector; 

1 9 determining an upper limit correction value for each symbol, the upper limit 

20 correction value being selectable to control an amount of signal to noise 

21 ratio deterioration; 
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22 applying a correction to the symbols to be transmitted in a direction opposite to 

23 the direction of the gradient vector of a magnitude not exceeding the 

24 determined correction values for each symbol; and 

25 transmitting the corrected symbols to the channel. 

1 2. The method of claim 1 wherein determining an upper limit correction value for each 

2 symbol, the upper limit correction value further comprises: 

3 computing an interpoint distance between symbols; 

4 selecting a correction value for a symbol as a value less than the interpoint 
p 5 distance to ensure that the symbol is not mistaken for other symbols. 

□ 1 3 . The method of claim 1 wherein applying a differentiable penalty function to the samples 

;-; w !2 having magnitudes exceeding the predetermined threshold comprises: 

; r 3 applying the function: 

b 4 

:| \{4k]-T) 2m if x[k] > T 

|| h(x[k]) = l 0 if \x[k]\<r 

f (x[k] + T) 2m if x[k]<-T 

6 where m is a positive integer that decides the severity of penalty, T is the 

7 predetermined threshold, x is the frame of data symbols expressed by: X = 

8 (r 0 , r x exp^), r 2 expG'6 2 ), . . expOe N/2 ^),r N/2 ) ? where r { and Q { denote 

9 the magnitude and phase of symbol in channel i, and k is the number of 
10 the symbol. 

1 4. The method of claim 3 wherein the net penalty function comprises: 
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N-l 



f(x) = "£h{x[k]) 



k=0 



1 5 . The method of claim 4, wherein the gradient vector is computed as: 



df _ £3 rfft(*[A:]) 

9r/ *=o 



cos 



+ 0. ;ze{l,...,tf/2-l} 



3f = g ^(xw) . a/ = g ^(xW) coc( _ ±) 

d£ S dx[k] U V 



1 6. The method of claim 1 wherein the gradient vector is computed only as a function of the 

□2 magnitude of the sample values. 

1=3 1 7. The method of claim 1 wherein computing a net penalty function value comprises adding 

^ 2 together the individual penalty function values computed for the samples having magnitudes 

MS! 

3 exceeding the predetermined threshold to generate the net penalty function value. 

fU 1 8. In a multicarrier communication system in which a signal to be transmitted comprises 

?D 2 data bits to be converted into a symbol modulated by each subcarrier of the signal prior to 

3 transmission on a channel, for a signal having a single peak in a frame, a method for minimizing 

4 a peak to average power ratio while minimizing introduction of errors into the signal to be 

5 transmitted: 

6 sampling the symbols to be transmitted of the frame; 

7 compare magnitudes of the samples of the frame to a predetermined threshold to 

8 determine whether sample magnitudes in the frame violate the 
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9 predetermined threshold, the predetermined threshold being selectable to 

1 0 control the number of samples violating the threshold; 

1 1 determining an upper limit correction value for each symbol, the upper limit 

12 correction value being selectable to control an amount of signal to noise 

13 ratio deterioration; 

14 computing a peak reduction kernel responsive to the upper limit correction values; 

15 responsive to determining a sample magnitude does violate the predetermined 

16 threshold, applying the peak reduction kernel to the sample to reduce the 

1 7 peak of the frame; and 

■fi 8 transmitting the modified symbol. 

It™ 

]Il 9. The method of claim 8 wherein computing a peak reduction kernel responsive to the 

p 2 upper limit correction values comprises: 

□ 3 determining a phase component and an amplitude component of the upper limit 

[U 4 correction values; and 

^ 5 setting the phase component of the upper limit correction values to zero to ensure 

6 that the peak reduction kernel has its peak value at the first sample of the 

7 frame. 

1 10. The method of claim 8 wherein applying the peak reduction kernel to the sample to 

2 reduce the peak of the sample comprises: 

3 rotating the peak kernel by amount to ensure a peak of the peak reduction kernel 

4 coincides with a peak of the frame; 
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5 determining whether the peak of the kernel has a sign equal to a sign of the peak 

6 of the frame; 

7 responsive to the signs of the peaks of the kernel and the frame being equal, 

8 multiplying the peak of the frame by minus one; and 

9 adding the peak of the kernel to the samples to reduce the peak of the frame. 

1 11. The method of claim 7 in a system in which more than one peak may be present per 

2 frame, comprising the steps of: 

3 responsive to determining that a sample magnitude exceeds the predetermined 

q 4 threshold, applying the peak kernel to the sample, wherein the peak kernel 

= P 5 applied for each sample has a magnitude scaled relative to an extent the 

I ;I! 6 sample magnitude exceeds the predetermined threshold. 

7~ 1 12. The method of claim 1 1 wherein the scaling factors are chosen to ensure a sum of the 

2 magnitudes of the kernels applied is equal to one. 

ft! 
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ABSTRACT OF THE DISCLOSURE 

A gradient sensitivity based method is described for reducing Peak to Average power 
Ratio in multicarrier communication systems. The method involves definition of a differentiable 
penalty function which has non-zero values whenever there is a peak violation. The gradient of 
the function with respect to the symbol magnitudes and phases is evaluated at the given symbol 
vector. Symbol perturbations are computed in the opposite direction of the gradient in such a way 
that the peaks are reduced and the amount of errors introduced in the symbols is not large enough 
to affect the decoding process. In order to reduce computation, a variant of the method is 
proposed where a peak reduction kernel is precomputed and stored. The kernel is rotated to the 
location of the peak and added with proper sign to reduce the peak to average power ratio. 
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